void copy_one(double *aout, double *ain, int n)
//Copies ain[1..n] to aout[1..n].
{
  int i;
  for (i=1;i<=n;i++)
    aout[i]=ain[i];
}
